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ABSTRACT 


The  content  of  an  ab-initio  time-dependent  theory  of  quantum  molecular  dynamics  of 
electrons  and  atomic  nuclei  is  presented.  Employing  the  time-depei^nt  variational  prindple 
and  a  family  of  approximate  state  vectors  yields  a  set  of  dynamical  equations  iq)pFQximating  tiw 
time-dependent  SduOdinger  equatitm.  These  equaticms  govern  the  time  evolution  of  the  relevant 
state  vector  parameters  as  molecular  orbital  coefficients,  nuclear  positions  and  momenta.  This 
q>proach  does  not  impose  the  Bom-Oppenheimer  q>proximati{»i,  does  not  use  potential  raergy 
surfaces  and  takes  into  account  elecmm-nuclear  coupling.  Basic  ccmsovation  laws  are  Bdly 
obeyed.  The  simplest  model  of  the  thetny  employs  a  single  determinantal  state  for  tiie  electrons 
and  classical  nuclei  and  is  implemented  in  the  cmnputer  code  ENDyne.  Results  from  this  ab- 
initio  theory  are  reported  for  icm-atom  and  ion-molecule  coUisitms. 
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Introduction 

Time-dependem  theoretical  methods  are  emer^g  important  tools  for  the  study  of  molecular 
processes.  Explicit  consideration  of  the  time  parameter  makes  it  possible  to  fdlow,  say,  the 
detailed  evolution  of  a  reactive  event  from  reactants  to  products  and  leads  naturally  to  the 
theoretical  determination  of  transition  probabilities  and  rates.  A  limiting  fact  for  most  dynamics 
methods,  both  time-dependent  and  time  independent  tmes  is  that  they  are  based  on  electnmic 
pcMential  energy  surfaces.  The  generatitm  of  a  quality  surface,  let  alone  several,  from  first 
principles  is  a  daunting  task  and  cmly  possible  for  small  molecular  systems.  One  often  resets 
to  the  construction  of  empirical  and  semi-empirical  surfaces  as  the  diatomics  in  molecule  (DIM) 
surfaces.  Although  inexpensive  to  generate,  such  surfaces  are  not  accurate  for  the  study  of 
detailed  dynamics.  If  only  a  single  such  surface  is  used,  as  is  often  tiie  case,  then  the  dynamics 
of  the  electrons  is  neglected.  Only  by  incmporating  two  or  more  surfaces  and  their  couplings 
is  the  effect  of  electronic  dynamics  included  [1]. 

There  are  i^yptoaches,  such  as  the  Qose-Coupling  method  [2,  3]  used  in  study  of  atomic 
and  molecular  collisicms,  where  electronic  dynamics  is  considered.  However,  these  methods 
onplcty  ncm-dynaittical  straight  line  or  Coulomb  nuclear  trajectnies  and  ignore  die  details  of 
tile  nuclear  dynamics.  While  fixed  trajectories  work,  well  for  collisicmal  phenomena  at  energies 
above  IkeV/amu,  they  lead  to  spurious  dynamics  at  lowm-  collisimi  energies. 

Another  body  of  work  has  its  origin  in  the  tqiproach  to  simulating  annealing  by  Car  and 
Pairinello  [4].  This  approach  has  been  generalized  to  study  dynamics  oi  molecular  systmns  or 
dusters  [S,  6,  7]  and  uses  a  fictitious  mass  associated  to  the  molecular  orbital  coefBcients  to 
generate  newtonian-like  equations  for  them  in  time.  The  method  attempts  to  follow  the  lowest 
pottmtial  energy  surftice  of  the  system,  without  having  to  perform  an  electnmic  optimizaticm  at 
each  new  nuclear  geometry.  Thus  it  is  very  much  like  dynamics  on  a  single  potmitial  energy 
surface,  but  without  the  cost  of  generating  tiie  full  surface. 

The  END  (Electron  Nuclear  Dynamics)  theory  of  Ohm,  Deumms,  et  al  [8, 9, 10, 11]  uses 
the  time-dqiendent  variational  prindple  Cn>VP)  to  generate  equaticms  that  qiproximate  tiie 
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tune-dependent  SchiOdinger  equation  and  govern  the  time  evolution  of  electronic  and  nuclear 
dynamical  variables  oa  a  generalized  phase  space.  The  choice  of  trial  wavefimction  for  the 
decttons  and  the  nuclei  determines  the  level  of  iq>proximati(m.  Both  electrcms  and  nuclei  ate 
treated  fully  dynanucally  including  coupling  terms.  The  dynamics  of  participating  nuclei  and 
electrons  is  subject  to  the  instantaneous  forces  due  to  the  full  Coulrunb  intnacticms  and  there  is 
no  need  to  generate  potential  energy  surfaces  arrd  associated  statitmary  states. 


A  related  method  by  Runge  and  Micha  [12,  13]  uses  the  Frenkel  variational  principle  to 
derive  equatitms  for  the  electtonic  degrees  of  fieedmn  and  tiien  use  the  eiktmal  a^rproximation 
for  the  nuclei  moving  in  die  potential  given  by  the  electnms  and  the  nuclear  Coulomb  repulsion. 
This  method  has  been  applied  to  proton-hydrogen  and  hydrogen-hydrogen  collisions. 

Ihe  simplest  level  of  the  END  theory  employs  a  sin^e  detnminant  for  the  electrons  and 
die  classical  limit  for  the  nuclei  This  level  of  qrproximation  has  been  implemented  in  a 
computer  code  ENDyne  [14].  In  spite  of  its  simplicity  this  level  of  theory  yields  accurate 
results  for  transition  probabilities,  differential  and  total  cross  sections  for  a  great  variety  of  ion- 
atom  and  ion-molecule  reactive  collisions.  The  END  theory  for  a  multiconfiguratitHial  and  for  an 
antisymmetrized  geminal  power  (AGP)  electronic  wavefimction  has  also  been  published  [15, 16]. 

The  next  section  gives  an  overview  of  die  END  theory  and  the  conservatitm  dieorems.  T 
Sectitm  in  discusses  the  effect  of  the  nonadiabatic  coupling  terms.  If  diese  couplings  are  ^ 
neglected,  the  calculation  tends  to  follow  a  potential  energy  surfisce,  in  a  Car-Psrinello-like 
fashitm.  An  example  shows  that  dynamics  without  these  terms  are  in  disagremnent  with  ~ 
rntperiment  When  die  couplings  are  included  agreement  is  exceUent 
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The  details  of  the  END  theory  [9]  and  the  associated  ENDyne  code  [11]  for  die  simplest  level 
of  theory  are  {Hiblished,  so  it  suffices  hett  to  give  a  summary  survey  to  khow  the  fiindamaital 
generality  and  pleasing  structure  of  the  dynamical  equations. 
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The  TDVP  [17]  for  a  full  quantum  i^tproach  requires  stadooarity  of  the  acticm  with  the 
Lasransian 

where  //  is  the  full  molecular  Hamiltonian  and  the  electronic-nuclear  many~body  wavefunction. 
The  derivative  ^  acts  cm  the  bra.  The  choice  of  trial  wavefunction  determines  the  details  of 
the  dynamical  equations.  Representing  the  nuclei  by  Gaussian  wavepackets  in  the  limit  of  zero 
widths  results  in  a  classical  description  with  nuclear  positicms  R  and  mcanenta  P  equivalent  to 
employing  the  Lagrangian 

^  ;  E  (ft  ® 

where  the  electronic  Hamiltonian  H^i  includes  the  nuclear-nuclear  repulsion.  The  END  dieoty 
even  admits  a  mixed  approach  with  some  of  die  nuclei  treated  by  classical  and  others  by  quantum 
mechanics. 


The  ^  quantum  TDVP  approach  with  a  completely  general  fonn  of  trial  funcdcm  yields 
the  time-dependent  Schrddinger  equation.  Restricting  the  trial  wavefuncdcm  to  a  particular  fonn 
leads  to  a  set  of  coupled  first-order  differential  equaticms  [17»  9]  governing  the  time  evolutkm  of 
die,  in  general,  complex  wavefunction  parameters.  The  parameters  assume  the  role  of  dynamical 
variables  and  fonn  a  generalized  phase  space.  This  phase  space  is  not  flat  but  has  a  nonunit 
metric  and  results  in  a  Hamilttmian  fonn  of  dynamical  equations.  The  chmce  of  wavefunction 
paranmters  bectmies  critical  to  ensure  nonredundancy  and  ctmtinuous  trajectories. 

Fot  the  case  of  a  single  determinantal  N-electrcm  state  vector  Tbouless  [18]  has  provided 
the  suitable  parametrization.  A  Thouless  determinant  \z)  is  an  antisymmetric  product  of 
iKmofthogonal  spin  orbitals 

K 

=  (3) 

a«N+i 

which  in  END  are  erqnessed  in  terms  of  an  attmiic  orbital  basis  set,  {^i;  i»l  JC],  centered  <m  the 
nudel  The  molecular  oiintal  coefficients  {z^}  and  dieir  complex  ocmjugates  are  tte  electronic 
dyiuinical  variables.  This  form  of  wavefunction  has  die  capacity  to  become  any  detmminantal 
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of  determinants  at  geometries  R  and  R!. 


The  END  theory  yields  a  dynamics  thtu  satisfies  the  conservation  laws  of  important  physical 
quantities,  as  total  energy,  total  momentum,  and  total  angular  momentum  [9, 1 1].  The  important 
conservation  of  total  linear  and  angular  mmnentum  only  holds  when  the  electron-nuclear  coupling 
terms  are  included. 


Electron-Nuclear  Coupling 

The  dynamical  metric  of  Eq.  (5)  couples  the  electronic  and  nuclear  degrees  of  freedom. 
The  thud  row  of  the  metric  matrix  has  cocq;>ling  terms  that  ensure  the  conservation  of  total 
mtxnentum  [9].  Tte  theory  must  satisfy  siKh  conservation  laws  to  be  able  to  produce  accurate 
detailed  information  about  nuclear  and  electronic  quantities.  One  might  think  that  at  higher 
energies  the  momentum  of  the  nuclei  overwhelms  the  contribution  from  the  electrons  when 
pr(q>erties  are  computed,  but  that  this  is  not  the  case  is  shown  by  our  results. 

The  first  row  of  the  metric  matrix  couples  the  electronic  degrees  of  freedom  to  the  nuclear 
velocities.  These  coupling  terms  make  up  tire  matrix  representation  of  tire  translatioD  qrerator  and 
provide  tire  transformation  of  the  z-coefficients  as  the  centers  of  the  basis  fiincticms  move  from 
(xre  instant  to  the  next  One  could,  alternatively,  use  a  basis  with  explicit  electron  translation 
factors  [20]  exp(-tvx  *  0*  where  tu  is  tire  velocity  of  center  but  in  its  full  implementatitm 
this  would  require  extensimi  of  current  integral  codes. 

Time-independent  approaches  introduce  ntmadiabatic  coupling  by  expanding  the  full  wave- 
fiinction  in  terms  of  adiabatic  states  with  each  term  being  a  product  of  an  mgenstate  of  the 
elecmmic  Hamiltonian  and  a  nuclear  wavefunctitm.  The  coupling  terms  arise  fitnn  tire  nuclear 
kinetic  energy  tqrerator  acting  on  tire  electronic  adiabatic  states,  which  are  parametrically  de¬ 
pendent  on  the  nuclear  positions.  Such  coupling  terms  are  absent  when  tire  full  wavefimction 
is  er^anded  in  terms  of  diabatic  states  witii  each  term  being  a  product  of  an  elecmmic  state 
independmit  of  nuclear  coordinates  and  a  nuclear  wavefunction.  Similarly,  in  time-dqrendent 
treatments  coupling  terms  emerge  from  using  elecmmic  basis  functitms  that  dqrend  oa  nuclear 
coordinates,  and  possibly  also  on  nuclear  velocities.  These  coupling  terms  are  absent  whmi  tire 
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etectronic  basis  is  fixed  in  space,  as  is  the  case  for  numerical  grids  fixed  in  the  lab  frame.  The 
dynamical  metric  matrix  in  Eq.  (S)  would  in  such  a  case  have  no  coupling  terms. 

The  influence  of  the  nonadiabatic  coupling  terms  in  Eq.  (5)  is  in  many  cases  essential  in 
order  to  describe  even  a  qualitatively  correct  dynamics  [21].  Neglecting  the  coupling  toms  in 
END  generates  dynamical  equations  identical  to  those  used  by  others  [22].  Other  iq)i»Daches 
that  neglea  etectron-nuclear  coupling  inclutte  dynamics  on  a  single  potential  energy  surface  and 
dynamics  using  fictitious  kinetic  energy  terms  to  force  the  dynamics  close  to  a  surface  as  does 
the  method  of  Car  and  Parrinello  [4,  6,  7].  Part  of  die  coupling  effects  can  be  handled  via 
electron  translation  factors,  as  is  done  by  Micha  and  Runge  [13],  but  that  still  accounts  for  (xily 
part  of  the  coupling.  A  dynamic  treatment  such  as  the  Qose  Coupling  metiiod  [23,  2]  uses 
electrcm  translation  factors,  but  limits  the  dynamics  to  fixed  trajectories,  which  also  neglects  full 
electron  nuclear  coupling  terms. 

Total  cross  sectitms  may  not  be  too  sensitive  to  the  neglect  oS  electron  nuclear  coupling 
terms,  but  differratial  cross  sectitms  and  details  such  as  the  presmrce  or  absence  of  rainbow 
scattering,  state  to  state  transition  probabilities,  branching  of  reactitm  channels,  etc.  cannot  be 
predicted  without  accounting  for  the  full  coupling. 

In  Fig.  1  tire  election  transfer  probability  at  250  eV  for  fee  H  system  as  a  function  of 
scattering  angle  is  depicted.  Results  from  two  experiments  [24,  25]  are  compared  to  the  END 
using  a  pVDZ  basis.  Experiments  have  a  finite  angular  resolutkm.  In  Ref.  [24]  die  resolutkm 
is  reported  to  be  Ad«0.4®  while  Ref.  [25]  reports  END  calculations  were  averaged 

to  reflect  the  finite  angular  resolutkms  of  fee  two  experiments.  Hiis  was  dcme  by  calculating 
probabilities  sufficiendy  close  together  (at  some  points  every  0.02**  ipait)  and  then  averaging 
die  probabilities  inside  ±  Ad/2  of  a  given  angle  $.  It  is  interesting  to  iKite  that  die  experimental 
values  reported  in  Ref.  [24]  are  actually  more  consistent  with  an  angular  resolution  of  1” 
rather  dum  0.4".  The  results  by  Helbig  and  Everhart  are  in  excellent  agreement  wife  END  for 
die  maxima,  but  the  minimum  value  is  too  high.  This  is  most  likely  due  to  some  nqierimental 
problnn  since  die  oscillations  should  be  symmetric  about  a  probalnlityttf  03  [24].  Other  sources 
of  error  in  the  eiqierimental  results  are  the  detmmination  of  die  dissociatitm  fraction  (H/Ha  in  the 
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odlision  region)  and  the  calibraticMi  necessary  between  the  different  pixKesses  (elastic,  charge 
tiansfer).  Another  problem  is  diat  as  the  scattering  angle  grows,  the  number  of  scattered  particles 
falls  very  stBq>ly.  This  actually  shifts  the  minima  and  maxima  towanls  larger  angles,  as  seen  in 
die  results  for  Htniver  et  al  at  angles  larger  than  4".  This  is  because  more  particles  are  measured 
omning  from  the  than  range  abmit  a  given  angle  0.  With  these  consitterations,  one 
can  otmclude  that  the  agreement  between  theory  and  experiment  is  excellent 

In  Fig.  2  the  END  results  for  the  same  transfer  probability  (using  perfea  resolutkm)  is 
ctHnpaied  to  the  results  that  completely  negle^  die  electron  nuclear  coupling  and  tme  for  which 
the  couplings  in  the  electronic  part  only  (t.e.  the  first  row  of  die  metric  of  Eq.  (S)  is  retained. 
The  latter  means  that  there  is  no  conservatitm  of  total  momentum.  It  is  clear  that  the  neglect 
of  the  coupling  terms  leads  to  wrcxig  behavim-  with  increasing  scattering  angle;  either  incorrect 
number  of  oscillations  or  shifted  positions  of  the  peaks. 

Reactive  Collisions 

ENDyne  calculations  of  total  and  differential  cross  sections  of  ion-atom  collisions  such  as 
protcm/hydrogen,  protcm/helium,  proton/lidiium,  alpha/helium  [9, 26, 27, 28]  widi  rather  modest 
Gaussian  basis  sets  have  shown  excellent  agreements  with  an  abundance  of  eiqierimental  data. 

For  instance,  proton/helium  collisions  have  been  studied  [26]  in  great  detail  for  a  variety 
of  collmon  energies  from  5(X)  to  5,000  eV.  The  classical  differential  cross  sectitms  for  electron 
transfix’  and  elastic  processes  are  calculated  using  a  pVDZ  basis  set  with  excellent  agreement  with 
experiment  for  rainbow  angles.  The  differential  elastic  cross  sections  with  a  simple  semiclassical 
correction  are  indistinguishable  from  the  experimental  results. 

Ongoing  work  for  small  collisional  systems  at  energies  from  a  few  to  a  fraction  of  an 
eV  yields  similarly  encouraging  results.  Since  END  imposes  no  restrictions  on  the  dynamics, 
beycmd  dua  of  basis  sets  and  wavefrinction  form,  basically  all  channels  accessible  from  given 
initial  conditions  can  be  studied.  The  analysis  of  processes  involving  icmization  with  (me  or 
more  dectton  departing  are  just  being  started  Processes  involving  km-molecule  coUisicms  with 
bond  breaking  and  bond  formaticm  are  readily  accessible  with  END. 
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PrcMon-hydiogeii  molecule  collisions  at  30eV  have  been  the  subjea  of  detailed  experimental 
invesdgatioa  [29]  .  Thme  have  been  several  attempts  to  explain  the  experiments  u^g  semi- 
classical  [29,  30]  and  quantum  descriptions  of  the  nuclei  [31,  32],  In  all  these  ^tfuoaches  DIM 
surfiuies  were  used  Florescu  et  al  modified  the  surfaces  to  achieve  quantitative  agreement  with 
experiment  [30].  Preliminaxy  results  from  END  show  that  the  rainbow  angle  is  calculated  more 
accurately  than  by  using  DIM  surfaces  and  the  i^proximations  inherent  in  the  'Hrajectory  Surface 
Htqyping  (TSH)  [29]  and  infinite  order  sudden  approximation  (lOSA)  [31,  32]  iq>proacbes. 

Experiment  puts  the  primary  rainbow  angle  for  charge  transfer  in  the  p+H2  collision  at  30eV 
between  6**  and  V.  The  TSH  and  quantal  lOSA  approach  methods  place  it  at  10°  [29,  31. 
32].  Preliminary  END  results  place  the  primary  rainbow  to  be  near  8.5°,  with  new  calculations 
lowering  this  result  The  primary  rainbow  angle  for  elastic  and  vibrationally  excited  collisions 
is  also  observed  to  be  between  6°  and  7°.  The  other  approaches  show  a  shift  of  approximately 
2J°  to  larger  angles,  while  preliminary  END  results  show  a  value  of  7.5°  for  the  rainbow  angle. 

Table  1  shows  the  different  channels  for  the  proton  on  H2  collision  as  identified  by  ENDyne 
(no  differentiation  among  different  vibrational  states  is  made  here).  It  is  clear  from  these  results 

Ihble  1  Product  chaimeb  for  various  oricntatioas  and  their  impact  porameta  ranges  for  IT  v-  H2  coUisioos  at 
30eV.  6  is  the  angle  formed  between  the  iiutial  velocity  vector  of  the  proton  (defined  as 
the  z  axis)  and  the  Ha  axis.  Ibe  initial  orientation  of  the  Ha  molecule  with  respect  to 
die  X'Z  plane  is  Ibe  impaa  parameter  b  is  zero  at  the  CM  of  Ha  and  increases  along 


the  X  axis.  A  range  of  probabilities  is  shown  for  the  charge  transfer  channeL 


(e.4i) 

Dissociation 

Proton  exchange 

Quurge  transfer 

(probabili^) 

(0.0) 

- 

0.0  <  b<  0.3 

0.3  <  b  (0.4-0.0) 

(45.0) 

0.65  <  b  < 

0.0  <  b<  0.6 

0.0  <  b<  0.65 

1.2  <  b  (0.24-0.0) 

(90.0) 

1.1  <  b<  1.35 

0.6  <  b<  1.1 

1.4  <  b  (0.26-0.0) 

(90.45) 

0.0  <  b<  0.4 

* 

0.4<b(0.2-0.0) 

that  the  assumptions  made  in  the  lOSA  and  the  use  of  a  DIM  surface  for  this  system  are  somewhat 
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questionable.  Tbe  orientation  oi  tiie  system  is  important  in  defining  the  outcome  of  a  collision. 
It  is  more  important  for  small  impact  parameters  titan  for  large  ones.  Furthermore,  dissociatimi 
is  ntX  accounted  for  in  the  lOSA  approach. 

The  total  breakup  channels  are,  of  course,  not  accessible  by  the  close-coupling  theory  or 
any  other  approach  using  the  so  called  impulse  approximation  that  does  not  permit  the  nuclear 
arrangement  of  the  molecule  to  adjust  during  the  collision.  Some  charmels  for  the  smallest 
impact  parameters  correspond  to  nuclear  transfer  and  should,  strictly  speaking  not  be  included  in 
the  electron  transfer  cross  section.  It  is  almost  impossible,  however,  to  separate  experimentally 
these  processes  from  the  transfer,  since  the  outgoing  proton  or  hydrogen  attnn  have  almost  the 
same  energy  as  the  incoming  proton. 

Therefore  it  is  essential  to  have  a  fully  dynamical  theory  accounting  for  electron  nuclear 
coupling  as  the  END  in  order  to  get  all  accessible  product  channels  and  to  interpret  the  experiment 
properly.  Artificial  restrictions  on  the  electron  nuclear  dynamics  by  prescribed  trajectories,  by 
frozen  nuclear  configurations,  by  forcing  the  dyrtamics  to  take  place  on  a  potential  energy  surface, 
or  on  several  surfaces  with  ad  hoc  hopping  introduce  nonphysical  dynamics  and  errors  that  are 
hard  to  ccmtrol. 

Another  system  under  study  is  the  low  energy  collisions  of  H2'*^  with  H2.  This  system 
exhibits  much  of  the  complexity  of  a  general  molecular  collisional  system  and  there  are  no 
accurate  potential  energy  surfaces  available.  For  the  few  theoretical  dynamics  studies  that  have 
been  performed  on  this  system  [33,  34]  approximate  potential  energy  surfaces  are  used,  and  the 
agreement  with  e]q)eriment  is  wanting.  Preliminary  results  for  END  using  a  Is,  2s  basis  for 
collinear  initial  conditions  with  a  single  impact  parameter  (Le.  assuming  ^herically  symmetric 
molecular  species)  show  good  agreement  with  experimental  total  cross  section  for  the  formatiem 
of  H3'*'  R  Further  studies  are  being  carried  out  beyond  these  simple  assumptions  for  the  initial 
conditions  including  all  channels. 
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Figure  CapUoQS 


Figtne  1  'Dinsfer  protebUiQr  v«mu  tcaaering  angk  for  p4-H  ooUision  m.  2S0  eV.  END  otiDg 
pVDZ  basis  and  three  avenged  results  using  an  angular  resolutiGn  ±A9  given 
in  parenthesis.  Experiments:  Expl:  Ref.  [24],  Exp2:  Ref  [25]. 

Figure  2  Transfer  piobabtlity  versus  scattering  angle  for  pt-H  ctdlision  at  250  eV.  usiiiiS  pVDZ 
basis.  Solid  line  same  as  in  Fig  1.  Dotted  line  without  any  electrQn>nuclear 
ooiq)ling.  Dashed  line:  n^lecting  momentum  conserving  couplings. 
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